Accepted to ApJL, 2011 Oct 25 

Preprint typeset using I^'T^]X style emulatcapj v. 11/10/09 



METALS REMOVED BY OUTFLOWS FROM MILKY WAY DWARF SPHEROIDAL GALAXIES* 

Evan N. Kirby^'^, Crystal L. Martin'', and Kristian Finlator^"' 
Accepted to ApJL, 2011 Oct 25 

ABSTRACT 

The stars in the dwarf spheroidal satelhte galaxies (dSphs) of the Milky Way are significantly more 
metal-poor than would be expected from a closed box model of chemical evolution. Gas outflows 
likely carried away most of the metals produced by the dSphs. Based on previous Keck/DEIMOS 
observations and models, we calculate the mass in Mg, Si, Ca, and Fe expelled from each of eight 
dSphs. Essentially, these masses are the differences between the observed amount of metals present in 
the dSphs' stars today and the inferred amount of metals produced by supernovae. We conclude that 
the dSphs lost 96% to > 99% of the metals their stars manufactured. We apply the observed mass 
function of Milky Way dSphs to the ejected mass function to determine that a single large dSph, like 
Fornax, lost more metals over 10 Gyr than all smaller dSphs combined. Therefore, small galaxies like 
dSphs are not significant contributors to the metal content of the intergalactic medium. Finally, we 
compare our ejected mass function to previous X-ray measurements of the metal content of the winds 
from the post-starburst dwarf irregular galaxy NGC 1569. Remarkably, the most recent starburst in 
that galaxy falls exactly on the ejected mass-stellar mass relation defined by the Milky Way dSphs. 
Subject headings: galaxies: dwarf — galaxies: abundances — galaxies: evolution — intergalactic 
medium — Local Group 



1. INTRODUCTION 

In a closed box model o f che mical evolution 
(jvan den Berghl 119621 : iSchmidtl I1963D . the peak of 
the stellar metallicity distribution lies at the nu- 
cleosynthetic yield. For most supernovae (SNe), 
the nucleo synthetic yield is abou t solar metallicity 
or hi g her (IWooslev fc Weaveii 119951 : ILimongi fc Chieffil 
120031: INomoto et all I2006D . However, the peaks of 
the stellar metallicity distributions of all known dwarf 
spheroidal galaxies (dSphs) are less than sola r (e.g., 
lAaronson fc Mouldl Il985l : iBuonanno et al.l 119851) . The 
dSphs likely lost metals during star formati on, which 
reduced their effective nucleos ynthetic yields (jHartwickl 
119761: lYoshii fc Arimotol 119871) . Almost no dSph with 
250 kpc of the Milky Way (MW) contains any de- 
tectable gas at present (e.g., iBlitz fc Robishawl 120001 : 
iGrcevich fc PutmanI 12009?) . Therefore, the missing met- 
als are not hiding in the galaxies' interstellar medium. 
The metals have been lost. 

Supernov a winds drive gas out of galaxies, particularly 
small ones (iLarsonl 119741: IDekel fc SiigiT986l) . It is pos- 
sible to detect this gas flowing out of galaxies, but it 
is very difficult to measure the metallicity or even the 
mass of galactic winds. Absorption line studies of distant 
galaxies, including galaxies at cosmological reds hift, can 
measure column densities and outflow vel ocities (" Martin 
[200llRuDke^raII[2005tlWeiner et al.l[200 9: Stcidc l et al. 
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[Mot iRubin et al.l IMol) but not mass or composition. 
Smaller galaxies may be even more prone to gas outflows 
due to their shallow gravitational potential wells. IMartinI 
(|1999l ) showed that the observed hot gas around gas-rich 
dwarf irregular galaxies (dirrs) in the nearby universe 
must be escaping from those galaxies in the form of out- 
flows. Just a few such galaxies are suitable for X-ray 
spectroscopy to measur e the compositio n of the escaping 
gas (e.g., NGC 1569 .Martin et al.ll2Ciol . 

Dwarf galaxies, which are numerous and have shallow 
gravitational potential wells, may be important contribu- 
tors to the metal content of the the intergalactic m edium 
(|Oppenheimer fc Davel [2008t iMartin et all [201^) or to 
their host galaxies. However, most nearby dwarf galax- 
ies, especially dSphs, have not experienced appreciable 
star formation for several Gyr. As a result, it is even 
more difficult to observe gas expelled from Local Group 
dwarfs because it must have happened so long ago. 

Kirby_et al. (2010, 2011b) recently measured the stel- 
lar metal content of dSphs. The amount of metals pro- 
duced over cosmic time for an individual dSph may be 
estimated from its stellar mass, an assumed stellar initial 
mass function (IMF), and theoretical SN yields and ex- 
plosion rates. The difference between the inferred metal 
production by SNe and the current stellar metal content 
gives the amount of metals the galaxy has lost. This 
letter presents that calculation and discusses the role of 
dSphs in enriching the intergalactic medium (IGM) and 
their host galaxies. 

2. CALCULATIONS OF METALS EXPELLED 

IKirbv et al.l (|2011al hereafter Klla) constructed a ba- 
sic chemical evolution model in order to explain the 
metallicity and [a/Fe] distributions of eight MW dSphs: 
Fornax, Leo I, Leo II, Sculptor, Sextans, Draco, Canes 
Venatici I, and Ursa Minor. The model assumes that 
each galaxy begins its life with a certain amount of gas 
(Mg(0)) but no stars. Primordial, metal- free gas flows 
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into the galaxy at a prescribed inflow rate {dMg/dt — 
Aint e~*/'^'"). Stars form according to a star formation 
rate law {dM^/dt = A^M^). Type II SNe explode 
from massive (10-100 Mq) stars, and Type la SNe ex- 
plode according t o an empirical delay time distribution 
(jMaoz et al.ll2010l ). Most importantly for this letter, the 
galaxy lose s a fix ed amount of gas (Aout) during each SN 
explosion. IKllal varied the parameters Mg{0), Ain, A^, 
Aout, Tin, and a until the modeled metallicity and [a/Fe] 
distributions matched the observed distributions. 

This model provides good matches to the abundance 
distributions of the less luminous dSphs. The inferred 
star formation histories also were qualitatively consistent 
with histories measured from color-magnitude diagrams. 
Essentially, nearly all stars in the less luminous dSphs 
are ancient, and they formed in 2 Gyr or less. However, 
the models for the most luminous dSphs (Fornax, Leo I, 
and to some extent, Leo II) have chemically measured 
star formation durations of only ~ 1 Gyr. The color- 
magnitude diagrams for these galaxies rule out such short 
star formation lifetimes, and they also rule out single 
bursts of star formation. These dSphs experienced sev- 
eral epochs of an increase in star fo rmation rate (SFR) 
followed by a decrease. The IKllal models permit only 
one episode of star formation (a single rise and fall of 
SFR). Therefore, the models are imperfect descriptions 
of Fornax, Leo I, and possibly Leo II. 

On the other hand, for an old enough population, such 
as the MW dSphs, most of the SNe that will ever ex- 
plode have already exploded. Therefore, the star for- 
mation history (SFR as a function of time) has a small 
effect on the total metal production integrated over cos- 
mic time. Because SN yields are metallicity-dependent, 
the star formation history does m atter t o a small degree. 

As a supplement to the work of lKllai we compute the 
mass in a number of metals. Mejected) that the dSphs 
lost to SN winds. The IKllal chemical evolution model 
tracks the amount and composition of expelled gas with 
each 1 Myr time step. The model was terminated when 
the dSphs were exhausted of gas. Although dSphs lose a 
small fraction of stars to tidal interactions with the MW 
(jMaiewski et al.|[2000l ). we neglect this deduction from 
the final stellar ma ss. The existence of the lumino sity- 
metallicity relation (|Mateolll998l:IKirbv et aLll2011b ^ lim- 
its the severity of tidal strip ping of luminous matter (also 
see iPenarrubia et"al]|2008() because tidally stripped ver- 
sions of massive satellites would be just as metal-rich as 
other massive satellites, but small, metal-rich satellites 
do not exist. 

In order to mitigate the uncertainties imposed on 
-^ejected from Uncertain star formation histories, we ad- 
vance the chemical evolution model for each dSph to an 
age of 10 Gyr. We compute the amount of ejecta — 
mostly in Fe — from Type la SNe over 10 Gyr, even 
though star forrn ation en ds in the mod el s long before 
then. Like IKllal we use INomoto ct al.l's (120061) nucle- 
osynthetic yields for Type II S Ne. Iwamoto et al.l's ()1999f ) 
yields for Type la SNe, and iMaoz et aLP s H2O1O0 delay 
time distribution for Type la SNe. The first Type la 
SN explodes 100 Myr after the birth of its progenitor, 
and the Type la SN rate declines as t~^'^. Based on 
the theoretical yields and IMF (Kr oupa et al.llT993l ) we 
have adopted. Type la SNe produce 94% of the iron over 



10 Gyr. Asymptotic giant branch stars are included in 
the chemical evolution model, but at these low metallic- 
ities, they are responsible for only a few percent of t he 
total production of Mg, Si, Ca, and Fe (Karakas 2 010f) . 

Table [1] gives A/cjcctod for Mg, Si, Ca, and Fe for 
the eight MW dSphs in our sample. We do not at- 
tempt to estimate errors on Mejected because the dom- 
inant source of uncertainty is the S N yields. Differences 
betw e en different yield m odels (e.g.. iWooslev fc Weaveii 
IT995t INomoto et al.ll20M) can be a factor of several for 
some elements. 

If Type la SNe have a prompt component 
(jScannapieco fc BildstenI |2005[ ). then the Type la 
SNe at early times could feed the same wind as the 
Type II SNe. However, the delayed Type la SNe at 
late times explode at a low rate. Without continuous 
explosions, it may not be possible to drive a wind. 
Nonetheless, nearly all dSphs are free of gas. There- 
fore, the "missing metals" — the difference between the 
current metal content of the dSph and the amount of 
metals estimated to have been produced by all of the 
dSph's SNe — have left the dSphs, whether in coherent 
w inds o r otherwise. 

IKllal assumed that the composition of the gas lost to 
SN winds was the same as the gas within the galaxy 
at the time of loss. However, outflows could be metal- 
enhanced, especially because the SNe that expel gas 
from the galaxy are the same ev e nts that enrich the 
galaxy with metals (jVadeil 119861 : iMac Low fc Ferraral 
[l999.) . The mass of the out flow is related inversely to 
its unknown metallicity, and IKllal showed that increas- 
ing the outflowing metallicity decreased the total amount 
of galactic mass loss required by as much as a factor of 
40. Given this uncertainty, we estimate only the mass in 
ejected metals and not the total mass of ejected material. 
The uncertainty in the metallicity of the outflows does 
not extend to the independent estimate of the ejected 
metal mass. 

3. RELATIONS WITH STELLAR MASS 
3.1. Total Metal Output 

Figure [1] shows the mass in each of several metals lost 
to SN winds summed over 10 Gyr. Not surprisingly, the 
dSphs with larger stellar masses eject more metals than 
the dSphs with smaller stellar masses. The Mejected — -M, 
relation is linear, which is the result expected if the dSphs 
lost all of their metals. The intercepts of the best-fit lines 
indicate that dSphs lost a constant 0.5% of their final 
stellar masses in the form of Fe, 0.2% as Si, 0.1% as Mg, 
and 0.01% as Ca. 

Of course, the dSphs could not have truly lost all of 
their metals because the stars are not metal-free. How- 
ever, the amount of metals retained (Mretained, the filled 
points in Figure [Ij in stars is tiny compared to the 
amount of metals produced. iKirbv et al.l (|2011br ) mea- 
sured the average metallicity of Fornax to be ([Fe/H]) — 
— 1.00 Given its stellar mass of 1.9 x 10^ Mq, Fornax's 
stars contain 4 x 10^ Mq of Fe. Even though Fornax 
retained the highest fraction of Fe of all eight dSphs con- 
sidered here, the galaxy presently has only 4% of the Fe 

^ We adopt the solar a bundance of Fc as 12 + log(n{Fe)/n(H)) = 
7.52 llSneden et aI|[T99g). We adopt 7.58, 7.55, a nd 6.36 for Mg, 
Si, and Ca, respectively lAnders fc Grevesselll989D . 
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TABLE 1 

Inferred Ejected Masses Over 10 Gyr, 



DSph 


TIT \ 

L (Lq) 


M* (Mq) 


MMg (Mq) 


Msi (Mq) 


n T f n T \ 
AiCa (Mq) 


Mpc (Mq) 


Fornax 


(1.8 ±0.5) X 10^ 


(1.9 ±0.5) X lO'^ 


1.2 X 10* 


2.9 X 10* 


1.9 X 10-^ 


8.3 X 10* 


Leo I 


(5.6 ± 1.7) X 10*^ 


(4.6 ± 1.4) X 1013 


3.7 X 10^ 


7.9 X 103 


4.9 X 102 


2.1 X 10* 


Sculptor 


(2.3 ± 1.1) X 10*^ 


(1.2 ±0.6) X 10*^ 


9.9 X 10^ 


2.1 X 10^ 


1.3 X 10^ 


5.5 X 10^ 


Leo II 


(6.6 ± 1.9) X 10" 


(1.4 ±0.4) X 10*^ 


1.2 X 10^ 


2.7 X 10^ 


1.7 X 10^ 


7.1 X 10^ 


Sextans 


(4.1 ± 1.2) X 10^ 


(8.5 ±2.4) X 10^ 


6.6 X 10^ 


1.5 X 10^ 


9.2 X 10^ 


3.8 X 10^ 


Draco 


(2.7 ±0.4) X IQS 


(9.1 ± 1.4) X 10^ 


7.0 X 10^ 


1.6 X 10^ 


9.7 X 10^ 


4.0 X 10^ 


Can. Ven. I 


(2.3 ±0.4) X IQS 


(6.3 ± 1.1) X 10^ 


4.9 X 10^ 


1.1 X 10^ 


6.8 X 10^ 


2.8 X 103 


Ursa Minor 


(2.2 ±0.7) X 10^ 


(5.6 ± 1.7) X 10^ 


3.9 X 10^ 


9.1 X 10^ 


5.5 X 10^ 


2.3 X 10^ 



References. — L (luminosity): Martin et al. (2008) for Canes Venatici L llrwin fc~H atzidimitriou ( 1995) 
otherwise. Af» (stellar mass): IWoo et al., (.2008 ). except that we assumed that Canes Venatici I the same 
M*/L as Ursa Minor. 



that it produced. In Ursa Minor, that fraction is only 
0.2%. 

The bottom panel of Figure[l]expresses metal retention 
in another way. The black circles show the stellar mass- 
iron m etallicity relation for MW dSphs (Kirbv et al] 
l2011bl) . The other points show the mass-metallicity re- 
lations for other elements. The smaller dSphs retain a 
larger fraction of a elements relative to iron than the 
larger dSphs. The dependence on stellar mass indicates 
that the larger dSphs formed stars for longer and there- 
fore experienced a higher ratio of Type la to Type II 
SNe. 

3.2. Relative Metal Output 

The mass-metallicity relation indicates that the more 
massive dSphs retained more of their metals in stars 
rather than losing them to galactic winds. Therefore, 
the less massive dSphs contribute a larger fraction of 
their metal mass to the IGM or to their host galaxies. 
To further explore the importance of small dSphs, we 
convolve the ejected mass function with the stell ar mass 
function for MW dSphs. IKoposov etahl (|2OO80 calcu- 
lated the luminosity function of MW dSphs by correcting 
the observed Sloan Digital Sky Survey luminosity func- 
tion for sky coverage and bias against low surface bright- 
ness dSphs. The resuh was dN/dMy = 10 x io01(a^v+5)_ 
Assuming the stellar mass-to-light ratio for an old pop- 
ulation, M^/L = 2, the relation becomes dN / dlogM^, = 
284 M""-^^ where is measured in Mq. 

Figure [2] shows the product of this dSph mass func- 
tion and the ejected mass function from Figure [1] Even 
though the small dSphs are more numerous, they are not 
significant contributors to IGM metals. Even if the rela- 
tion for the smaller dSphs is extrapolated to zero mass, 
then all dSphs the size of Leo I or smaller contributed 
just 7.5 X 10'' Mq of Fe to the IGM over 10 Gyr. Fornax 
alone provided more Fe, and even more massive galaxies 
likely make Fornax's contribution negligible. The same 
argument applies to the other metals as well. 

3.3. [a/Fe] 

Averaged over a long enough time, the ratio of the 
numbers of Type II to Type la SNe asymptotes to the 
same value in all dSphs as long as the IMF is the same. 
As a result, the [a/Fe] ratios of ejecta from dSphs in- 
tegrated over 10 Gyr are not functions of stellar mass. 
Expressed relative to the solar ratios, the 10 Gyr ra- 



tios for our adopted SN yields and our adopted IMF are 
[Mg/Fc] = -0.45, [Si/Fe] = -0.15, and [Ca/Fe] = -0.30. 

4. COMPARISON TO NGC 1569 

iMaver et al.l ()2006l ) theorized that gas-rich dirrs can 
be converted into gas-free dSphs by interaction with 
their host galaxies. Tidal stirring can turn disks into 
spheroids, and ram pressure stripping — in addition to 
supernova-driven outfiow — can remove all of the gas. 
With this relationship between dIrrs and dSphs in 
mind, we compare our inferences for the metals ejected 
from dSphs to the observed metal content of the wind 
presently flowing out the post-sta rburst dirr NGC 1569 
( M, = 1.8 X 1() 8 lh7^ [T988h . 

iMartin et al.l (|2002[ ) estimated that the most recent 
starburst in that galaxy ejected Mq ~ 3.4 x 10^ Mq 
of oxygen. They based this estimate on an X-ray spec- 
tral fit to several a elements, including O and Mg. In 
order to compare to our own ejected metal estimates, 
we note that Mq is related to A^Mg by a multiplicative 
factor (about 15, which is also the solar ratio) because 
both elements are formed in similar nucleosynthetic pro- 
cesses in massive stars. Therefore, the starburst must 
have ejected Mjvig ~ 2.3 x 10'^ M© of magnesium. The 
mass of stars formed in the burst was M» « 3.4 x 10^ Mq 
over 10-20 Myr. 

NGC 1569 is more massive than any of the dSphs we 
consider, and the SFR in the burst was much higher than 
the average rates in the dSphs. Nonetheless, the outfiow- 
ing magnesium mass from the recent burst in NGC 1569 
agrees remarkably well with the relation between M-^g 
and defined by MW dSphs. Figure [T] shows that the 
NGC 1569 outflow falls exactly in line between Leo I and 
Leo II. Of course, we must keep in mind that we are con- 
sidering only the very recent starburst and not the entire 
galaxy, which is 50 times more massive. However, only 
the ejecta of the current starburst is directly measurable, 
and the measurement is uncertain by a factor of at least 
2 

IMartin et al.l (|2002f ) also measured the [a/Fe] ratio to 
be 0.3-0.6 (2-4 times solar). This value is larger than 
what we infer for the MW dSphs. However, the winds 
flowing out of NGC 1569 came from a very short burst — 
likely short enough that Type la SN ejecta is not present 
in the winds. Therefore, the [a/Fe] ratios merely reflect 
the fact that the star formation was intense in the recent 
burst in NGC 1569 but not in the dSphs. 
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Fig. 1. — Top: The mass in four different metals inferred to have 
been lost from (hollow points) or retained in (filled points) MW 
dSphs as a function of their present stellar mass. The hollow points 
show the "missing metals," which should have been produced by 
the SNe in these dSphs but are not present in the stars. The masses 
of ejected metals include the ejecta of Type la SNe that exploded 
up to 10 Gyr after the cessation of star formation. The dashed lines 
illustrate slopes of AfJ"", which is the expected relation if dSphs 
lost all of the metals they produced. The large star shows the Mg 
mass in the current outflow from NGC 1569 versus the stellar mass 
of the most recent starburst (Martin et al. 2002). Botto m: The 
stellar mass-metallicity relation for the stars in dSphs (Kir by et al.l 
I2011bl) . Metallicity is given separately as [Fe/H], [Mg/H], [Si/H], 
and [Ca/H]. 



IMartin et al.l 's completely independent method for 
measuring the metal content of galactic winds has placed 
a different, though possibly related, type of dwarf galaxy 
on the same Afojcctod".^* relation as the one defined by 
MW dSphs. This result has several implications. First, 
NGC 1569 retained almost none of t he metals pro duced 
in its recent starburst, a result that IMartin et al.l found 
from their own chemical evolution model. Second, the 
theoretical SN yields are accurate enough to place the 
observed magnesium mass ejected from NGC 1569 on 
the same relation as the the masses ejected from dSphs, 
which were calculated from the theoretical yields. Fi- 
nally, even a galaxy as massive as NGC 1569 {a^ = 
21 km s-\ iStil fc Israeli 120021 compared to 10 km s'^ 
for a typical dSph) loses nearly all of the metals it pro- 
duces. 
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Fig. 2. — The metals ejected from dSphs convolved with the 
luminosity function for MW dSphs {dN/d\ogAh oc M~°-^^, 
IKoposov eral|[200a) . The dashed line shows AfO'^^ -pj^g 

axes are 

logarithmic, and the metal output from small dSphs is negligible 
compared to large dSphs. 

5. DISCUSSION 

The possibility that the majority of metals that cur- 
rently reside i n the IGM were originally ejected by low- 
mass systems (iQppenheimer fc Davel[2008t IMartin et al.l 
2010) emphasizes the need for constraints on models 
for galactic outflows. Galaxy evolution models (e.g., 
iSomerville et al.ll2001l : IQppenheimer fc Davell20M ) gen- 
erally assume that star-forming galaxies drive outflows 
whose metallicity tracks the interstellar medium, and 
whose outflowing mass is proportional to M^/tr". Here, 
fj is the velocity dispersion that scales with the halo mass 

1/3 

{a oc ), and a is a scaling parameter. Setting a = 1 
yields a momentum-driven scaling, and a = 2 yields a 
energy-driven scaling. Some models (e.g., iBower et al.l 
2006) invoke a stronger scaling (a — 3) in order to match 
the faint end of the luminosity function. Let us assume 
that the metallicity of outflowing gas Zw tracks the 
stellar metallici ty Z^,, which in tur n scales with stellar 
mass as M"-^^ (|Kirbv et al.l 1201 IbD . Let us further as- 
sume the stellar mass-halo mas s scaling that is observe d 
at low masses (M* oc Mg-^^ iLeauthaud et al.l 120111) . 
The mass in ejected metals is the outflow mass multi- 
plied by the wind metallicity. Therefore, the model as- 
sumptions correspond to the assumption that Mojoctcd oc 
iM,/a'^)Zw = Afi-si/M^"/' = ^^i.3i-a/6.5i^ ^^i.ig 
and M^ ^" for momentum-driven and energy-driven scal- 
ings, respectively. Taken at face value, our observation 
of the ejected mass relation, Mojoctcd oc M^-°°, favors an 
energy-driven scaling. 

We emphasize that these inferences depend strongly 
on the assumed stellar mass-halo mass relation as well 
as the metallicities within the outflows. They also ne- 
glect many physical processes, such as tidal stripping, 
that could preferentially affect low-mass systems. Fur- 
thermore, much of the metals from Type la SNe are 



Outflows from dSphs 



5 



lost wefl after the cessation of star formation and well 
after the galaxies are devoid of interstellar gas. That 
makes the distinction between energy- and momentum- 
driven winds poorly defined. These uncertainties em- 
phasize the need for improved observational constraints 
on the kinematics of low-mass systems as well as the 
mass-loading and metal mass fractions of outfiows from 
low-mass galaxies. 

The biggest limitation of our approach at measuring 
.^ejected IS that wc havc measured the mass only in met- 
als. We provide no constraints on the total outflowing 
mass, including hydrogen. Therefore, we also provide no 
constraints on the metal fraction of the outflows. One 
way to constrain the metallicity is to fix the total mass 
participating in gas flows and star formation. For in- 
stance, we could assume that the baryonic mass associ- 
ated with a dSph obeys the cosmic baryon-to-dark mat- 
ter ratio. Unfortunately, t he disconnect between stel- 
lar a nd dark matter mass (iMateol 119981: iStrigari et al.l 
|2008[) and the tidal interactions (jPeharrubia et al.ll2008lf 
of MW dSphs make it impossible to estimate their total, 
original dark matter masses. 

Semi-analytic models (e.g.. iCooper et aP l201Clf ) may 
provide a path forward. By treating a variety of pro- 



cesses that likely affected dSph growth, such as tidal 
stripping, pre-enriched inflows, and returning outflows, 
they provide a flexible but complete account of how 
outflows impact metallicities. A model galaxy whose 
chemical abundance s match observed abundance distri- 
butions (jKirbv et al. 2011a,b) and our inferred ejected 
metal masses might inform us about these processes. We 
encourage interested modelers to interrogate their exist- 
ing semi-analytic models for unique constraints that our 
observations place on the way in which low-mass galax- 
ies process their gas into outflows, stars, and evolving gas 
reservoirs. 
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